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Abstract 

The global regularity for the incompressible magnetohydrodynamic equations (MHD) 
in three dimensions is a long standing open problem of fluid dynamics and PDE theory. 
The Navier-Stokes equations can be viewed as a special case of MHD with a constant 
magnetic field, whose global regularity problem is known as a Clay Millennium Prize 
Problem. In this article, we prove the global regularity of axially symmetric solutions to 
the ideal MHD in three dimensions for a family of non-trivial magnetic fields. The proofs 
are based on the special structures of MHD and can of course also applied to the resistive 
MHD. Our result might indicate that there are richer fantastic research topics in MHD 
than Navier-Stokes equations. 

Keywords: Magnetohydrodynamics, global regularity, axial symmetry. 

1 Introduction 

Magnetohydrodynamics (MHD) is to study the behavior of an electrically-conducting fluids. 
Examples of such fluids include plasmas, liquid metals, salt water, etc. The field of MHD 
was initiated by Hannes Alfven, for which he received the Nobel Prize in Physics in 1970. 
However, the mathematical theory on MHD is still very little known until today. 

The fundamental concept behind MHD is that magnetic fields can induce currents in 
a moving conductive fluid, which in turn creates forces on the fluid and also changes the 
magnetic field itself. MHD owes its peculiar interest and difficulty to this interaction between 
the field and the fluid motion. The set of equations which describe MHD are a combination of 
the Navier-Stokes equations of fluid dynamics and Maxwell's equations of electromagnetism. 

We focus on the incompressible case. Let B denote the magnetic field, E the electric field, 
u the bulk plasma velocity, J the current density, p the plasma pressure and [i > the fluid 
viscosity. The momentum equation is 

d t u + u • Vu + Vp = fiAu + J x B. (1.1) 

The ideal Ohm's law for a plasma is given by 

E + u x B = 0. 
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Faraday's law is 

d t B = -V x E. 

In an imperfectly conducting fluid the magnetic field can generally move through the fluid 
following a diffusion law with the resistivity of the plasma serving as a diffusion constant. 
This includes an extra term in Ohm's Law and gives the resistive MHD. The Faraday's law 
in this case reads 

dtB = V x (u x B) + uAB. (1.2) 

Here v > is the resistivity. The low-frequency Ampere's law neglects displacement current 
and is given by 

J = —V x B. (1.3) 
Mo 

Here hq > is the magnetic constant. The incompressibility and the magnetic divergence 
constraints are 

V-u = 0, V-B = 0. (1.4) 

The MHD is a combination of ([TT]) . (fL2|) . (fL3|) and ([Ti]) (see, for instance, [9]). 

In this paper, we are going to make a preliminary study on the global well-posedness 
of general solutions to MHD equations. Our main result concerns the following three- 
dimensional ideal MHD: 

d t u + (u • V)u + Vp = Au + (V x B) x B, 
< 3,B = Vx(ux B), (1.5) 
V-u = 0, V-B = 0. 

We have set the constant no and /i to be 1. For the ideal MHD, the fluid has so little resistivity 
that it can be treated as a perfect conductor. This is the limit of infinite magnetic Reynolds 
number. For applications of ideal MHD, see, for instance, 

We will show the global existence of classical solutions to the 3D incompressible ideal 
MHD in the axially symmetric case for a class of non-trivial magnetic fields. The theorem is 
stated as follows: 

Theorem 1.1. Suppose that uo and Bq are both axially symmetric divergence-free vectors 
with Uq = and Bq = Bq = 0. Moreover, we assume that (uo,Bo) E H s with s > 2 and 

— f 1 £ L°°. Then there exists a unique global solution (u, B) for the ideal MHD (11.5P with the 
initial data (uq,Bo) which satisfies 

t 5 

||u(t,.)||i 2 + ||B(t,.)||^+ / \\Vuf H2 ds<e eet \ 

Jo 

The notations used here can be found in section 2. Note that the Paraday's equation 
for B in (11. 5j) is exactly the same as the vorticity equation for the 3D incompressible Euler 
equations (by identifying B and V x u). This may lead to an essential difficulties for the 
global well-posedness of the ideal MHD (II. 5h in general case. Indeed, the global regularity of 
(jl.5p is even widely open in the two-dimensional case if the magnetic field is non-trivial. We 
achieve Theorem 11.51 by exploring the underlying special structures of the MHD system in 
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axially symmetric case. The magnetic stretching term B • Vu in Paraday's equation can be 
absorbed into the convection term by dividing the equation by r. This yields that II = — 
is only transported by the velocity field u. On the other hand, by dividing r in the vorticity 
equation, one can absorb the vortex stretching term into the convection term, leaving only 
one term involving II as a forcing one in f2 equation. 

The above result is motivated by a novelty observation on the following resistive MHD: 

' d t u + (u • V)u + Vp = Au + (V x B) x B, 

< 9 t B = AB+Vx(ux B), (1.6) 
V-u = 0, V-B = 0. 

v. 

Here we have set all the constants /xo, A 4 , v to be 1. We have the following theorem: 

Theorem 1.2. Suppose that uq and Bo are both axially symmetric divergence-free vectors 
with Uq = and Bq = Bq = 0. Moreover, we assume that (uo, Bo) € H and — E L°° . Then 
there exists a unique global solution (u, B) for the resistive MHD (jl.6p with the initial data 
(uo,Bo). Moreover, (u, B) is smooth in the sense that (u(i, •), B(t, •)) £ H s for any s > 
and t > 0. 

As is well-known that any smooth solution of 3D Navier-Stokes equations is automati- 
cally a solution of MHD. In particular, the solution to 3D axially symmetric Navier-Stokes 
equations without swirl (see [8j [IT]) also solves 3D MHD. However, we believe that there 
are richer structures in the case of MHD so that there are more examples of global solutions 
to MHD equation. For instance, u = B = e* A uo is a solution of (|1.6p with the associated 
pressure p = — ■j-- The above Theorem 11.61 gives another solution to (|1.6|) with axially 
symmetric initial data whose magnetic field is non-trivial. Our result might indicate that 
there are richer fantastic research topics in MHD than Navier-Stokes equations. We hope 
that this could initiate further study on this topics along this line of research. 

Our key observation is that if we rewrite the 3D incompressible axially symmetric Navier- 
Stokes equations in terms of u = u r e r +u z e z and b = u s eg, then there is only a sign differencial 
between the Navier-Stokes equations (|2.9h for (u, b) and the resistive MHD (|1.6j) for (u, B) 
(see Remark 12.11 in section 2 for details). As a consequence, there is an extra maximum 
principle for the quantity n of the resistive MHD (jl.6p . while for Navier-Stokes equations, 
the maximum principle only holds for the quantity T = ru e . This sign difference finally 
results in the global well-posedness in Theorem 11.61 for 3D MHD, while the corresponding 
problem for 3D Navier-Stokes equations is still widely open. 

At last, let us consider the perfect resistive MHD: 

'd t u + (u • V)u + Vp = (V x B) x B, 

< d t B = AB + V x (u x B), (1.7) 
V-u = 0, V-B = 0. 

Unfortunately, at present, we are unable to establish a global regularity result for fjl . 71) under 
similar settings as in Theorem 11.51 and Theorem 11.61 But we can prove the existence of weak 

1 In fact, the pressure is also changed. But the pressure is not a troublesome term for our purpose due to 
the divergence-free condition. 
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solutions. We remark that such a result in the general case is also widely open. Our theorem 
is: 

Theorem 1.3. Suppose that uo and Bo are both axially symmetric divergence-free vectors 
with Uq = and Bq = Bq = 0. Moreover, we assume that uo G H l n W 2 '2, B G H 1 and 

y ei 2 fl L°° . Then there exists a global weak solution (u, B) G L°°(0, oo; L 2 ) for the perfect 
resistive MHD fjl .6f) with the initial data (uo,Bo). 

Due to Theorem 11.21 it is expected that we can get a family of approximated solutions 
by adding the viscous term eAu into the momentum equations and then apply the viscosity 
vanishing method to yield the desired weak solutions to (jl.7p . The point is to get enough 
a priori estimate so that we can deal with the Lorentz forcing term when taking the limit 
e — ¥ 0. For a technical reason, we require that the initial data have relatively high regularity. 
It is possible to obtain the result under less regularity of initial data, but we don't pursue 
that in this paper. 

Before ending the introduction, let us mention some important results in the field of 
incompressible MHD. The local well-posedness of the resistive MHD (|1.6|) was established in 
|13j where the authors also proved the global well-posedness in 2D case. A nontrivial blowup 
criterion for the perfect resistive MHD fj LTD was established in terms of only L| (BMO) norm 
of vorticity of the velocity field in [10J . Recently, Lin, Xu and Zhang [11] made an important 
progress on the global well-posedness of classical solutions for the 2D ideal MHD (II. 5p under 
the assumption that the initial velocity field and the displacement of the magnetic field from 
a non-zero constant is sufficiently small in appropriate Sobolev spaces. Cao and Wu [2] 
proved the global regularity of 2D resistive MHD with partial viscosity and resistivity. We 
also emphasis the partial regularity theory and Serrin type criterions in [H [6] , and various 
blowup criterions in [1] [3] (see also the reference therein) . 

The remaining of this paper is simply organized as follows: In section 2 we will derive the 
axisymmetric MHD in cylindrical coordinate. We will also make a comment on the difference 
between the resistive MHD fll .6D and the axially symmetric Navier-Stokes equations and 
prove a maximum principle for n. We will prove Theorem II .11 in section 3. In section 4 we 
present the proofs for Theorem 11.21 and Theorem 11.31 

Acknowledgement. The author would like to thank professor FangHua Lin, professor 
Horng-Tzer Yau and professor Shing-Tung Yau for their interest in this work. 



2 Axially Symmetric MHD and A Maximum Principle 

In this section we will first derive the incompressible axially symmetric MHD in cylindrical 
coordinate. Then we show that the quantity n satisfies a maximum principle. We also present 
an interesting connection between the axisymmetric MHD studied in Theorem 11.21 with the 
axisymmetric Navier-Stokes equations with non-trivial swirl u® (see (|2.6f) and (|2.9|) ). 

Let us begin with some notations. A point in M 3 is denoted by x = (x\, X2, z). Let 
r = yjx\ + x\ and 





















(i: 








2) 






v 5 I 







4 



be the three orthogonal unit vectors along the radial, the angular, and the axial directions 
respectively. An axially symmetric solution to the 3D incompressible MHD (jl.6p is a solution 
(u, B,p) which takes the following form 

f ~ 

u(t, x) = u r (t, r, z)e r + u (t, r, z)eg + u z (t, r, z)e z , 
B(t,x) = B r (t,r,z)e r + B 9 (t,r,z)e e + B z (t,r, z)e z , 
^p(t,x) =p(t,r,z). 

We will also write the vorticity field V x u in cylindrical coordinate: 

V x u(i,x) = u> r (t,r, z)e r + oj (t,r,z)eg + u z (t,r, z)e z , 

where 



Define 



u r = -d z u , J = d z u r -d r u z , co z = -d r (ru ). 



n = — , n = — , v = ru 6 



By expanding the Lorentz force term as 



(V x B) x B = B • VB 



IBI 



(2.1) 



and then taking the inner product of u and B equations with e r , e$ and e z , respectively, we 
can derive the resistive MHD in cylindrical coordinate: 

d t u r + u r d r u r + u z d z u r - ^ + d r P 

= (A-^)u r + B r d r B r + B z 8 z B r - 

d t u e + u r d r u e + u z d z u e + ^f- 

= (A - + + B z d z B e + S^L, 

dtu z + u r d r u z + u z d z u z + d z P 

= Au z + B r d r B z + B z d z B z , 



d t B r + u r d r B r + u z d z B T 



(A - \)B r + B r a r u r + 5 2 &u r , 



d t B e + u r d r B e + B z d z u e + 2^1 

= (A - + B r <9 r ?/ + B*0*tt* + 

fl^S* + u r <9 r .B 2 + u z d z B z 

= AB Z + B r d r u z + £ 2 <9 z u 2 , 



where the pressure is given by 



P = p + 



IBI 



(2.2) 



(2.3) 
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The incompressible constraints are 



v r B r 
d r u r + — + d z u z = 0, d r B r + — + d z B z = 0. 



(2.4) 



The general axially symmetric resistive MHD is governed by (|2.2p and (|2.4p . In this paper, 
we consider a family of solutions with the form 



u(i,x) = u r (t, r, z)e r + u z (t, r, z)e z , B(i, x) = B (t, r, z)e$. 



(2.5) 



It is easy to check that (u 9 , B r , B z ) can be zero for all time if they are zero initially. In this 
case, (u, B, P) in (|2.5|) and (|2.3|) is governed by 



d t u r + u r d r u r + u z d z u r + d r P = (A - 



r l ) r 



d t u z + u r d r u z + u z d z u z + <9 2 -P = Aii z , 
«9 tJ B e + n^ r 5 e + u^ 2 5 e = (A - ^)B 9 + 

together with the incompressible constraint 

V 

d r u r + — + cw = 0. 
r 

To avoid the explicit presence of pressure, we also need the vorticity formula of (|2.6p 
d t B 9 + u r d r B 9 + u z d z B 9 = (A - \)B 9 + 



(2.6) 



(2.7) 



d t u 9 + u T d r u 9 + u z d z oj 9 + £ 



A -4, W 



iv,( &(b 9 ) 2 



(2.8) 



Remark 2.1. It is well-known that the axially symmetric Navier-Stokes equations (in the case 
of B = 0) are (see, for instance, 



r 



<9f-u r + u r d r u r + u z d z u r + d r p = (A — ^-)u r + 
c^-u 2 + u r d r u z + u z d z u z + c\p = Au 



+ u r <9 r ?/ + u z <9 2 -u 9 = (A 



If we denote u = u r e r + u 2 e 2 and b = u 9 gq, we can rewrite the above axially symmetric 
Navier-Stokes equations as 



d t u + u • Vu + Vp = Au - b • Vb, 
d t b + u • Vb = Au - b • Vu, 
V-u = V-b = 0. 



(2.9) 



If we compare the MHD equations (jl.6p with the Navier-Stokes equations (|2.9p . we find that 
we can recover (II. 6p from ()2.9p by changing the sign of the terms b • Vb and b • Vu. The 
significance of the tiny difference, especially, the sign of b • Vu, yields a much stronger a 
priori estimate in the MHD case. 
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Proposition 2.1 (Maximum Principle). Assume that (u, B,P) is a smooth solution to (12. 6h 
with or without resistivity and r,z)\ — > as r + \z\ — > oo /or eac/i /ixec? i > 0. Then the 
quantity U satisfies the maximum principle 

l|n(t,-)IU» < ||n(o,-)IU», vt>o. 

Proof. Let us first consider the case of zero resistivity. By dividing the equation for B e by r, 
one has 

d t Il + u r d r n + u z d z U = 0. (2.10) 
Since all functions are axially symmetric, we can write fl2. lQf) as 

d t u + u • vn = 0, 

which gives the maximum principle for II in the case of zero resistivity. 
Similarly, in the resistive case, we have 

d t U + u r d r U + u z d z U = (d 2 + -d r + d 2 z )U. (2.11) 

r 

Clearly, for smooth solutions, one has H(t,0,z) = for all time. Consequently, | IX | can not 
reach its maximum on the axis r = if it is not identically zero. Suppose that |IT(io , to, zo)| = 
su Po<t<t Oll^ 00 a t some to > 0- Then one must have ro > 0. Now one can restrict the 
equation (|2.11|) on the domain 

Q e = {(t,x):0<t<t , e<r < -,\z - z \ < -}. 

By viewing (|2.1ip as a five-dimensional draft-diffusion equation, we conclude by strong max- 
imum principle that II = ||n(0, Olli 00 on Qe- Taking e — > 0, one has II = 1 1 TT(0 , Oll-k 00 for all 
< t < tQ. Since H(t, 0, z) = 0, one in fact has II(t, •) = for all time, which contradicts to 
the fact that II doesn't vanish identically. Hence, one has \\T(t, Oll-L 00 < ||r(0, Olli 00 fo r an Y 
t > 0. 

□ 



3 Proof of Theorem 11.11 

In this section we prove Theorem ll.il Throughout this paper, we will use A\ < A2 to denote 
that A\ < C0A2 and A\ ~ A2 to denote that Cq 1 A2 < A\ < C0A2 for a generic positive 
constant Co > 1 and two positive quantities A\ and A2- 

Proof of Theorem \l.l\ Let us rewrite the vorticity equation in (|2.8p in terms of 

d t n + u r d r n + u z d z n = (M + -d r + d 2 z )n - d z u 2 . 

r 

By taking the L 2 inner product of the above equation with f2 and preforming the standard 
energy estimate, one has 

1 d 



2dt 



\\n\\ 2 L2 - J n(d 2 + ^d r + d 2 )ndx 
^ J ( u r d r n 2 + u z d z n 2 )dx - J no z u 2 dx. 
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Using the incompressibility condition (12. 7p and the fact of dx = 27rrdrdz, one has 

J (u r d r n 2 + u z d z Q 2 )dx = 



and 



z)\ 2 dz. 



- I n(d 2 + -d r + d 2 z )ndx=\\Vn\\ 2 L2 +27r [ |O(i,0, 

J r JR 

By integration by part and interpolation, we have 

| J nd z u 2 dx\ < \\u\\ 2 L4 \\d z n\\ L2 < ^linilialinnl. + ^\\d z n\\h- 



Consequently, one has 



^||o||| 2 + ||v^||| 2 < iiniilaiiniUcc. (3.1) 



Similarly, using equation (|2.1ip and preforming the L 2 energy estimate, one has 

\\U(t,-)\\ L 2<\\U \\ L 2, Vt>0. (3.2) 
Consequently, by Proposition 12.11 and using (|3.ip . (|3.2p . we have 

Mt,-)\\ L 2 <i + Vt, f ||vn||| 2 di<i + t, vt>o. (3.3) 

J o 

Here we used that f2o G L 2 which is due to the fact that uo E H 2 and 

|V(V x u)| 2 = \ (e r d r + -e e d e + e z d z )oj e e e \ 2 = \Vuj 8 \ 2 + M 2 . 

T 1 

Similarly, one also has IIo £ L 2 since Bo G H 1 which is also used. 

To proceed, we need a technical lemma regarding the property of a Riesz operator on M 3 . 
We first recall the following weighted C alder on- Zygmund inequality for a singular integral 
operator with a weight function which is in the A p class (see Stein |14| pp. 194-217 for 
details). Let fC be a Riesz operator in W 1 and w(x) be a weight in the A p class (see page 
194 of |14j for definition). One can extend the C alder on- Zygmund inequality for the singular 
integral operator with the integral having weight function w(x). Specifically, for 1 < p < oo, 
there holds 

\\lCf\\iA>(jt») < ll/llip(Rn), V / G L p (K. n ). 
The following lemma plays an essential role in our global regularity analysis. 
Lemma 3.1. There holds 



T ! j-T 



i 



r-V(t,-)||r~dt< sup ||n(v)H! 2 / \\vn(t,-)\\l 2 dt. 

o 0<t<T Jo 

Remark 3.2. We pointed out that in [7] the authors have established an inequality ||?'~ 1 9 2 ii r '||iP < 
for 1 < p < oo for M 2 x T , where T 1 is a one-dimensional torus. 
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Proof. We follow the proof in [TJ. By the incompressible constraint (I2.7p . we can introduce 
the angular stream function ijj such that 



and 



r 



(3.4) 



We divide by r in (|3.4p . which gives that 



(3.5) 



Following [7], we interpret the Laplace operator in (|3.5p as a five-dimensional one. We 
formally write 



y = (yi,y2,V3,yi,z), r = yyj + y^ + yl + yj, A y = (d? + -d r + d%). 

This way we have %- = (-Ay)- 1 ^. In the remaining part of the proof of this lemma, we 
will use a subscript y to denote the derivatives with respect to y. 
It is clear that 



V ; 



(e r d r + -e 9 d e + e z d z )(e r d r — ) + V<9 2 — ® e 2 



,0 



1 n V>' 



e r g> e r d r h e e (g) e e -<9 r h (e 2 g> e r + e r <g> e 2 )<9 2r h e 2 <g e z <9: 



Consequently, one has 



' r r 



(3.6) 



On the other hand, one also has 



V- 6 



i> 6 



Vy — = (e r d r + V<? + e z d z )(e r d r — ) + V y d z — (g> e 



,0 



e r (g e r <9„ h Vee r d r h (e 2 (g e r + e r ® e 2 )<9 2r h e 2 <g> e 2 <9 

e r (g) e r <9 2 h (J - e r (g e r )-<9 r — 

/y* ly* ly* 

+ (e 2 (g e r + e r <g e 2 )5 2r h e 2 <g> e 2 <9 z — . 



where Ii 



o 



14x4 




and Vq is defined by 



Ve = V - e r (e r • V y ) - e z d z , e r = - 





fyi\ 




/ON 


1 


2/2 









2/3 


,e 2 = 





r 


2/4 









w 




w 
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Clearly, e r (g> e r , Iq — e r e r , e z ® e r , e r (g) e z and e 2 (g) are all mutually orthogonal. 
Consequently, one also has 

iv^i 2 -i<i 2 + iMi 2+ i<i 2+ i<i 2 <"> 

By (|3.6p and (|3.7p . we have 

iv^p, a r r^i'+ii^i'+i^i'+i^n'nw, 

i r i J-ooJo x r ' r r r 1 / 

= r r^i'+ii^i'+i^+i^nw** 

J-ooJo V r r r r r / 



-oo JO 



V 2 , — | P w(r)r 3 dr<i2; 



~ J \V 2 y (-A y )- l ^\ p w(r)dy, 

where w(r) is a weight function io(r) = -y. 

Let 1 < p < oo. Using Lemma 2 in [7] (see also a general version in Lemma 15.11 in 
Appendix of this paper), we have 

J \vl(-& y )-^ w ( r )dy < J \^\ P w(r)dy 



l^fdx. 
r 



\ V ^\P dx < f \^L\ p dx. (3.8) 



Consequently, one has 

Repeating the above procedure, one also has 

\V 2 ^\ P dx< fl^Fdx. (3.9) 



i 



Taking p = 2 in ([3SD and and using the interpolation inequality ||/||l«> < ||V/||| 2 ||V 2 /||£ 2 
in M 3 , one has 



r T r T 

/ Ur-VfcOlLootft = / \\r^d^ e (t r )\\ Loc dt 
Jo " JO 



< 



o 



' |V^(r-V e (t,-))||i||V 2 ^(r"V e (t,-))||! 2 ^ 
< sup ||Q(V))lll2 [ T \\dMt,-)\\hdt. 

0<t<T JO 

This finishes the proof of the lemma. □ 



10 



Now we derive an L°° estimate for B . Ignoring the viscosity in the equation of B in 
(12,61). one has 



\\B e (t,-)\\ L ~ < \\B$\\ Lao + f \\B e (s, .)|| £ oo II — 1| da. 

Jo r 

By GronwalPs inequality and using (|3.3p and Lemma 13.11 we have 

||£ e (V)||L~ < \\Bl\\ La °eti\\ r ~ ld * ur{s ^ LOods < e tl . (3.10) 
Let us coming back to (|2.8p and estimate that 

v r 

< ||y|| Loo / {uj ) 2 dx+\\Ii\\ L oo\\B e \\ L 2\\d z J\\ L 2 

< llylL- / (-*) 2 ^ + ^Piil-ii^ + ^ii^ni, 

Recalling the following basic energy law 

— (llulli, + ||B||| 2 ) + jT (||Vu||| 2 + ||VB||i 2 )d S = 0, (3.11) 
and using the a priori estimate in (|3.3p and Lemma 13. H one has 

[|Vxu(V)||l» <e tS , /V(Vxu)||| 2 ^<e 41 , V f > 0. (3.12) 

JO 

The next step is to bootstrap the regularity of u and B. We are going to show the 
L^fO, T], Lip(M 3 )) estimate of u. We will make use of the structure of the ideal MHD in (|2.6p 
to avoid some possible technical complications. The key observation is that we can write the 
vorticity equation as 

d t (V x u) + V • [(V x u) x u] = A(V x u) - d z (UB e e e ). 

Here by the maximum principle in Proposition ^. H and (|3.10p . one has T1B S £ L°°([0, t], L°°(IR 3 ) 
Moreover, we can apply (|3.12p to bootstrap the regularity of (V x u) x u. Then we may 
apply the standard parabolic estimate to get the L°°([0, t], L°°(R 3 )) estimate for V x u. 
We first perform L 4 energy estimate for (|4.4p and derive that 

T I* 

< \\ — \\ Lao {J)UX +\\U\\ L o \\B 9 \\ L o a \\d z \J\^ L 2\\J\\ L 2 

< II7 IL- / (^fdx + iliniii-HB^ili-iic/iii, + \\\d z \J\\\l2. 

Using the a priori estimate in (|3.3p and Lemma 13.11 one has 

III, 61 ! 2 !! 2 4- IIVL , e l 2 ll 2 < J 1 

II 1^ I Hz,°°([0,t], L 2 (M 3 )) II V l W I \\L^([0,t],L 2 (^ 3 ^ >' 6 ■ 
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By interpolation inequality, one has 

5 

1011 1011 t ^ 

\\ u llL°°([0,t],L 4 (M 3 )) + \\ u lli 4 ([0,t],L 12 (M 3 )) ~ e • 

On the other hand, by Sobolev imbedding, one also has 

2 

\\ u \\L^([0,t],L°°(R3)) < H u llL°°([0,t],L 2 (R 3 )) + ll w llL°°([0,t],L 4 (R 3 )) ~ e ■ 

Hence, we have 

5 

||(V x u) x u|| L 4 ([0it]iL i2 (IR 3 )) < e* 3 . 



Write 



V x u = e tA 



10 

A standard parabolic estimate gives that 



V x u - f e (*- s ) A (V • [(V x u) x u] + d z {UB e e e ))ds. 
Jo 



5 

|VV x u|| L 4 ( [ 0it ] ii i2( R 3)) < e . 



By Sobolev imbedding, we have 

5 

l|Vu|| i 4([o it ] ) Lo ( R 3)) < e . (3.13) 
Now let us derive the L 1 ([0, T], Lip(lR 3 )) estimate of B. We first write 



d t B + u • VB = — B. 

r 



Applying V, one has 



u r 1 
d t VB + u • VVB = -Vu • VB + — VB + Vu'Ue e + (V-> r B. 

Note that 

(V-)n r B = -—He r , 

one has 

||VB(t,-)IU« < ||VB ||l-+^ (||Vu||l<» + 11^11^) ||VB(a,-)|U»da 

ft u r 
+ / (||Vu|| X oc + — 1| )\\Il(s,-)\\ L °ods. 
Jo r 

We can use (13.3H . (13.13p . Lemma [37X1 and Gronwall's inequality to estimate that 



|VB(t,.)||L~ <e et '. (3.14) 
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The a priori estimates fj3. 13H and (|3.14j) are enough for the global regularity of the ideal 
MHD equations (jl.5p . Indeed, applying the standard H 2 energy estimate, one has 

^(l|V 2 u(V)||| 2 + ||V 2 B(t, .)||i 2 ) + ||V 3 u(V)|li 2 

< J (- V 2 uV 2 (u- Vu) + V 2 uV 2 (B • VB))dx 

+ J (- V 2 BV 2 (u- VB) + V 2 BV 2 (B- Vu))dx 

< ^||V 3 u(t, 0||! 2 + ||u(t, l|V 2 u(t, -)||i 2 + ||B(t, HV 2 B(i, .)||| 2 
+ (||Vu(t, + || VB(t, -)\\ L oo) (||V 2 B(t, -)||! 2 + ||V 2 u(t, OH!.)- 

Here we also used the Gagliardo-Nirenberg's inequality ||V/||^ 4 < ll/IU 00 ||V 2 /||^2, the inte- 
gration by parts and j V 2 B(u • V)V 2 Bcte = 0. Consequently, one has 

5 

\V 2 u(t,.)\\ L2 + \\V 2 B(t,.)\\ L2 <e eet \ 

V t > 0. 

f \WM\hdt<e^, 

We finished the proof of Theorem 11.51 □ 

4 Proofs of Theorem 11.21 and Theorem 11.31 

In this section we prove Theorem 11.21 and Theorem 11.31 
Proof of Theorem M.SX Similarly as in obtaining (|3. 1 j) . one has 

1 d f 

j^lMlia + HVnilia < | J Sld z Tl 2 dx\ 

< ||n|||^||np 10 ||^n|| i2 

LT 

< ||n|||^||n||f 2 ||vn|| L2 ||c^|| L 2. 



Consequently, one has 



d_ 

Jt 



^-\\n\\ 2 L2 + ||vo||| 2 < iiniil^iiniil.iivnii 2 ,. (4.1) 



Applying a similar argument to II equation in (|2.1ip . one has 

^||n||2 2 + ||vn||i 2 <o. (4.2) 

Clearly, the combination of (|4.ip . (|4.2p and the maximum estimate in Proposition 12.11 gives 
the following a priori estimate 



\\u(t,-)\\ L2 + \m,-)\\ L2 <i (Vi>o), / (||vn|| 2 2 + ||vo|| 2 2 )dt <i. (4. 

Jo 



3) 



13 



Now let us come back to the equation of u in (|2.8p . Applying the standard energy 
estimate, one has 



i|/i"T<fa+/(iv.r+^ 



(4.4) 



4(B") 2 ., 



u) dx. 



Using Sobolev imbedding theorem and interpolation, one has 



-dx 



< \\u r \\ L 2 ||n||i6 \\0J \\l3 



< iiuniL 2 |ivoii L2 |i^ii| 2 ||v^" 5 

i, 



liaii n L 2 

< IL.HI2 IIV70II2 , n, .0||2 



U 



L 2 



II 2 - 



On the other hand, it is clear that one also has 



I 



-u) dx 



< \\n.\\ LO o\\B°\\ I ?\\d g u a \\ ri 



L 2 + ^ll^z^ Hi 2 - 



Using the a priori estimate (|4.3p . the basic energy law ()3.11|) and Proposition 12. 1\ we have 
||Vxu(t,-)||L 2 <1 (Vi>0), / ||V(Vxu)||| 2 )dt<l. (4.5) 

JO 

The a priori estimate (|4.5p is enough to get the global regularity of the resistive MHD 
(|1.6p . Indeed, by using the equation of B, one can easily verifies that VB also satisfies (j4.5|) . 
We have finished the proof of Theorem 11.21 □ 

Next, let us present the proof of Theorem 11.31 

Proof. We first consider the following resistive MHD with e fluid diffusion which will serve 
as the regularized system of ([1.70 : 



d t u + (u • V)u + Vp = eAu + (V x B) x B, 
d t B = AB + V x (u x B), 
V • u = 0, V • B = 0. 

To impose an initial data for ([4.60 . we choose two non-negative mollifiers (pi £ Cq°( 
<fo G C^QR 2 ) such that 

J (f>x[z)dz = J (j>2{xi, X2)dx\dx2 = 1, 4>2 is symmetric. 

Denote x = (x±,X2,z), y = (yi,y2,z') and define 



(4.6) 



and 



Jef(x) 



</>l (e z)<t) 2 (e 2/1, e y 2 )f(x -y)dy. 
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For given initial data uo and Bo as in Theorem 11.31 define 

u = rWr- 1 *®]* + (J e u z )e z , B = r^r' 1 B e )]e e . (4.7) 

Clearly, due to our choice of J e , the regularized data in (|4.7|) is also axially symmetric. Since 
uo € H 2 , one has r _1 UQ = —d r UQ — 8 z Uq G iT . Moreover, since 

IV^Cr- 1 ^)]*)! 2 = IV^Cr-^JD^ + I^Cr-^)! 2 



+ |rV y ^(x-yJCr" 1 ^)^^]^! 2 
< I^Cr^l' + KrV^drl^KDl' + IVJeK 112 



01 > 



one has Uq G ff 1 for any fixed e > 0. Here we used the notation r p = \Jp\ + p\ for a point 
p € M 3 . Similarly, one also has Bg G -ff 1 for any fixed e > 0. It is obvious that ^P- G L°°. 
Hence, all the assumptions on the initial data in Theorem 11.21 are satisfied by u and Bo e . 
We can then apply Theorem [L2] to conclude that there exists a unique global smooth solution 
(u e ,B e ) to for each e > 0. 

We are going to take the limit e — > 0+. Note that for each fixed e > 0, by energy estimate, 
one has the following a priori estimate 



|(u%BO|| l2 <||(u ,B )|| l2 <||(uo,B )||hi, 
, |VBi L2a0jOo]iL2) < ||K,B )|| i2 < ||(uo,B )|| H i ; 

which are independent of e. Moreover, since 

1 d 



(4- 



and 



~l|n'& + ||vn<& = o, 



— 11^1% <-2 [ ^-^WdJWfdx < ||O e ||%|| n e|i 6 ||a 2 n e || L 2 , 
at / 



(4.9) 



one also has 

l ne ||L 2 + ll Vne |lL 2 ([0,oo),L 2 ) < ||n^|| L 2 < ll-BoH^i, 

mi L§ <l|fi ll zS + llvni|| 2([0)Oo))i2) <c, 

for a positive constant depending on ||uo||#2 and ||Bo||#i, but independent of e. Note that 
ft = = d z [J e (r- l u r Q )] - 9riJ ; uZ ° ] G W 2 '^ by Hardy inequality. 

For any R > and any vector <j) G Cq°(Br) with V • <j> = 0, taking V> G Co* 3 (-B3.fi:) with 
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ip = 1 on B2R, and using (14. 8p . (14,911 and (|4.6p . we can estimate that 



dtu e (bdx 



= -J x u e ) x u e K<pdx + J (V x B e ) x WR<j)dx 

- J (I - ^)(V x u e ) x u e K(f>dx 

= -J *prtt e e e x u € K<j)dx + J(V xB € )x B e K(j)dx 

+ J (1 - ^)u e (8) u e : VTl(j)dx + J u e ® u e : V(l - V) ® 
< 3i?||^|| L 3||(l-V)ul L 6||</»|| L6 + ||VB e || L2 ||Bl L6 ||«/»|| L 3 



+ Kill, I 



|x-y|>2i? 



|ac - y\ 4 l<A(y)|dy|| L2 + K 8) u £ || L2(£ . 2ii) ||c/>|| L 2 



IVB 



e II 2 



H 1 ■ 



Consequently, one can derive that 



Similarly, one also has 



\dtU e \\L]H-i(B R ) < C R- 



l^t^llijH-^Bji) ~ Ca- 



using the standard compactness theorem (see, for instance, [16]), one has 

(u e ,B £ )^(u,B) in L 3 ([0,oo),L? ocal ) 

for a sequence of e — > 0. By a diagonal argument, we get a subsequence of (u e ,B e ) such 
that the above strong convergence is true for any R > 0. Let r)(t, •) G C 1 ([0, oo), C^°) be 
an arbitrary test function with V • rj(t, x) = 0. We write the resistive MHD system in the 
following weak form 

'// u e (%7 + (u e • V)r/ + eAr/) - (B £ • Vr])Wdxdt = 0, 
/ / B e (d t rj + (u e • V)t/ + A??) - (B e • Vrj)u e dxdt = 0, 
JJ B e • Vr]dxdt = 0. 

Clearly, by taking the limit e — > 0+, one has 

'// u(d t r] + (u • V)r/ - (B • Vr])Bdxdt = 0, 
/ / B(d t r] + (u • V)t/ + At/) - (B • Vr])udxdt = 0, 
// B • Vr/dxcft = 0. 

Since r/ is divergence and arbitrary, one concludes that there exists a distribution p such that 
(u, B,p) satisfies (jl.7p in the sense of distribution with the initial given data and (u, B) 6 
L°°(0,oo;L 2 ), VB G L 2 (0,oo;L 2 ). 

□ 
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5 Appendix 

In this appendix we first prove that w(y) = r a is a A p for Riesz operator in M 5 under 
—4 < a < 4p(l — ~) . The case of a = — 2 has been studied in [7]. 

Lemma 5.1 („4 p Weight). Let 1 < p < oo and w(y) = r a , y£l 5 . Then w{x) is in A p class 
if-A<a< 4p(l - |). 

Proof. Recall that a real valued non-negative function w(x) is said to be in A p (W l ) class if 
it satisfies 

1 f 1 f ~ 

sup ^TgT y w(x)da^ (rgj y w(:e)~pc£e^ 9 < oo. 



Here p and g are conjugate indices with 1 < p < oo. 

For any ball -B C M 5 , denote B = B(yo,R). It is easy to see that if ro > 2i?, one has 
r ~ ro for any x £ B. Consequently, for any aGK, one has 



1 f 1 /* ^ 

— — / w{x)dx ) ( -j-=rr / pdx) q 
l-°l Jb ' M-°l ' 



5 (t4 /" rffdxVvi / /■„ ' //./■]" f 1. 



5i 



On the other hand, if ro < 2R, then for a + 3 > —1 and — ^ + 3 > —1, one has 







(w,L 










/TO 


<(- 






ro- 




ct _ 



1 /" , 4 _4 



B 



r +B /-3R „„ . _ . E 



r -R JO 



Noting that the condition on a is — 4 < a < 4p(l — |) , we in fact have completed the proof 
of the lemma. □ 

Next, let us provide a simpler proof of Theorem 11.11 under an extra condition on the initial 
data n G L 3+s . 

Theorem 5.2. Suppose that all the assumptions in Theorem \l.l\ are satisfied. Furthermore, 
if Qo G L 3+s for some 5 > 0, then there exists a unique global solution (u, B) for the ideal 
MHD (I1.6P with the initial data (uo,Bq) which satisfies 



\u(t, -)\\ 2 H2 + ||B(t, Ollfp + r ||Vu||^ 2 & < e< 

JO 



3 
.12 



Proof. We of course still have the a priori estimate (|3.3|) . On the other hand, by multiplying 
the equation of f2 by |f2| 1+<5 f2 and then preform the standard energy estimate as follows: 



|||0||^ + ||V|0|^||^ < I / n^8 z tfdx 



< ^l|v|o|^||| 2 + ||n||* A ||o||^ 

< ^liv|^|^|l! 2 + ||n||^||n||l+i||0|llt 5 . 
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Using (I3.3p . one has 

\m,-)\\ L3+ s <l + tK Vt>0. (5.1) 
Now let us recall a Biot-Savart law (see [15]): 

K(( , x)| </ ^„« + r[ ^ iy . 

J\y—x\<ir \ x U\ J\y—x\>Ar \ x V\ 



Here y = (yi,y2,z'). Denote r' = v2/i + v\- O ne nas 

l« (*>r,*)| < / ; | _ 2 dy + r ■ - dy. 

J\y—x\<4r r \ x U\ J\y—x\>Ar r \ x V\ 

Note that r' < \r — r'\ + r < \x — y\ + r. Hence, if \x — y\ < 4r, one has r' < 5r, and if 
|x — y\ > 4r, one has r' < 2\x — y\. Consequently, by Young's inequality, we have 

\u r (t,r,z)\ < r-7-rz * U < rWr-r^W 3 \\n\\ L s,i. 

\X\ A \X\ Z "L2< 

Here L p ' q denotes the usual Lorentz norm which is defined as 

(t\{x:\f(x)\>t}\- P ) q j) q , Zft" = 2£ = 

Noting that 

1 11 1 2 

Ti2 r*,- = supt|{x : — >f}|s < 1, 



and using the real interpolation and (13. 3p . (15. ip . we have 

\ U r (t,r,z) I 3(1 2 * 4) ,. .. 3^jy 1 

I " I ~ II"IIl 2 ll^llxa+a ^i + r 2 - 

Using (|5.'2p . the remaining proof can be carried out similarly as in section 3. □ 



(5.2) 
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